Search for: All records

We measure the thermal electron energization in 1D and 2D particle-in-cell simulations of quasi-perpendicular, low-beta (β_p= 0.25) collisionless ion–electron shocks with mass ratiom_i/m_e= 200, fast Mach number${\mathcal{M}}_{\mathrm{ms}}=1$–4, and upstream magnetic field angleθ_Bn= 55°–85° from the shock normal$\hat{\mathit{n}}$. It is known that shock electron heating is described by an ambipolar,B-parallel electric potential jump, Δϕ_∥, that scales roughly linearly with the electron temperature jump. Our simulations have$\mathrm{\Delta }{\varphi }_{\parallel }/(0.5m_{\mathrm{i}}{u_{\mathrm{sh}}}^2)\sim 0.1$–0.2 in units of the pre-shock ions’ bulk kinetic energy, in agreement with prior measurements and simulations. Different ways to measureϕ_∥, including the use of de Hoffmann–Teller frame fields, agree to tens-of-percent accuracy. Neglecting off-diagonal electron pressure tensor terms can lead to a systematic underestimate ofϕ_∥in our low-β_pshocks. We further focus on twoθ_Bn= 65° shocks: a${\mathcal{M}}_{\mathrm{s}}=4$(${\mathcal{M}}_{\mathrm{A}}=1.8$) case with a long, 30d_iprecursor of whistler waves along$\hat{\mathit{n}}$, and a${\mathcal{M}}_{\mathrm{s}}=7$(${\mathcal{M}}_{\mathrm{A}}=3.2$) case with a shorter, 5d_iprecursor of whistlers oblique to both$\hat{\mathit{n}}$andB;d_iis the ion skin depth. Within the precursors,ϕ_∥has a secular rise toward the shock along multiple whistler wavelengths and also has localized spikes within magnetic troughs. In a 1D simulation of the${\mathcal{M}}_{\mathrm{s}}=4$,θ_Bn= 65° case,ϕ_∥shows a weak dependence on the electron plasma-to-cyclotron frequency ratioω_pe/Ω_ce, andϕ_∥decreases by a factor of 2 asm_i/m_eis raised to the true proton–electron value of 1836.

 

We study a relativistic collisionless electron–positron shock propagating into an unmagnetized ambient medium using 2D particle-in-cell simulations of unprecedented duration and size. The shock generates intermittent magnetic structures of increasingly larger size as the simulation progresses. Toward the end of our simulation, at around 26,000 plasma times, the magnetic coherence scale approachesλ∼ 100 plasma skin depths, both ahead and behind the shock front. We anticipate a continued growth ofλbeyond the time span of our simulation, as long as the shock accelerates particles to increasingly higher energies. The post-shock field is concentrated in localized patches, which maintain a local magnetic energy fractionε_B∼ 0.1. Particles randomly sampling the downstream fields spend most of their time in low field regions (ε_B≪ 0.1) but emit a large fraction of the synchrotron power in the localized patches with strong fields (ε_B∼ 0.1). Our results have important implications for models of gamma-ray burst afterglows.

 

We perform non-radiative two-dimensional particle-in-cell simulations of magnetic reconnection for various strengths of the guide field (perpendicular to the reversing field), in magnetically dominated electron–positron plasmas. Magnetic reconnection under such conditions could operate in accretion disc coronae around black holes. There, it has been suggested that the transrelativistic bulk motions of reconnection plasmoids containing inverse-Compton-cooled electrons could Compton-upscatter soft photons to produce the observed non-thermal hard X-rays. Our simulations are performed for magnetizations 3 ≤ σ ≤ 40 (defined as the ratio of enthalpy density of the reversing field to plasma enthalpy density) and guide field strengths 0 ≤ Bg/B0 ≤ 1 (normalized to the reversing field strength B0). We find that the mean bulk energy of the reconnected plasma depends only weakly on the flow magnetization but strongly on the guide field strength – with Bg/B0 = 1 yielding a mean bulk energy twice smaller than Bg/B0 = 0. Similarly, the dispersion of bulk motions around the mean – a signature of stochasticity in the plasmoid chain’s motions – is weakly dependent on magnetization (for σ ≳ 10) but strongly dependent on the guide field strength – dropping by more than a factor of two from Bg/B0 = 0 to Bg/B0 = 1. In short, reconnection in strong guide fields (Bg/B0 ∼ 1) leads to slower and more ordered plasmoid bulk motions than its weak guide field (Bg/B0 ∼ 0) counterpart.

 

The recent discovery of astrophysical neutrinos from the Seyfert galaxy NGC 1068 suggests the presence of nonthermal protons within a compact “coronal” region close to the central black hole. The acceleration mechanism of these nonthermal protons remains elusive. We show that a large-scale magnetic reconnection layer, of the order of a few gravitational radii, may provide such a mechanism. In such a scenario, rough energy equipartition between magnetic fields, X-ray photons, and nonthermal protons is established in the reconnection region. Motivated by recent 3D particle-in-cell simulations of relativistic reconnection, we assume that the spectrum of accelerated protons is a broken power law, with the break energy being constrained by energy conservation (i.e., the energy density of accelerated protons is at most comparable to the magnetic energy density). The proton spectrum is${\mathit{dn}}_p/{\mathit{dE}}_p\propto E_p^{-1}$below the break and${\mathit{dn}}_p/{\mathit{dE}}_p\propto E_p^{-s}$above the break, with IceCube neutrino observations suggestings≃ 3. Protons above the break lose most of their energy within the reconnection layer via photohadronic collisions with the coronal X-rays, producing a neutrino signal in good agreement with the recent observations. Gamma rays injected in photohadronic collisions are cascaded to lower energies, sustaining the population of electron–positron pairs that makes the corona moderately Compton thick.

 

Magnetic reconnection is often invoked as a source of high-energy particles, and in relativistic astrophysical systems it is regarded as a prime candidate for powering fast and bright flares. We present a novel analytical model—supported and benchmarked with large-scale three-dimensional kinetic particle-in-cell simulations in electron–positron plasmas—that elucidates the physics governing the generation of power-law energy spectra in relativistic reconnection. Particles with Lorentz factorγ≳ 3σ(here,σis the magnetization) gain most of their energy in the inflow region, while meandering between the two sides of the reconnection layer. Their acceleration time is$t_{\mathrm{acc}}\sim \gamma {\eta }_{\mathrm{rec}}^{-1}{\omega }_{\mathrm{c}}^{-1}\simeq 20\gamma {\omega }_{\mathrm{c}}^{-1}$, whereη_rec≃ 0.06 is the inflow speed in units of the speed of light andω_c=eB₀/mcis the gyrofrequency in the upstream magnetic field. They leave the region of active energization aftert_esc, when they get captured by one of the outflowing flux ropes of reconnected plasma. We directly measuret_escin our simulations and find thatt_esc∼t_accforσ≳ few. This leads to a universal (i.e.,σ-independent) power-law spectrum${\mathit{dN}}_{\mathrm{free}}/d\gamma \propto {\gamma }^{-1}$for the particles undergoing active acceleration, and$\mathit{dN}/d\gamma \propto {\gamma }^{-2}$for the overall particle population. Our results help to shed light on the ubiquitous presence of power-law particle and photon spectra in astrophysical nonthermal sources.

 

The non-linear interaction between electromagnetic waves and plasmas attracts significant attention in astrophysics because it can affect the propagation of Fast Radio Bursts (FRBs) – luminous millisecond-duration pulses detected at radio frequency. The filamentation instability (FI) – a type of non-linear wave–plasma interaction – is considered to be dominant near FRB sources, and its non-linear development may also affect the inferred dispersion measure of FRBs. In this paper, we carry out fully kinetic particle-in-cell simulations of the FI in unmagnetized pair plasmas. Our simulations show that the FI generates transverse density filaments, and that the electromagnetic wave propagates in near vacuum between them, as in a waveguide. The density filaments keep merging until force balance between the wave ponderomotive force and the plasma pressure gradient is established. We estimate the merging time-scale and discuss the implications of filament merging for FRB observations.

 

Abstract In galaxy clusters, the intracluster medium (ICM) is expected to host a diffuse, long-lived, and invisible population of “fossil” cosmic-ray electrons (CRe) with 1–100 MeV energies. These CRe, if reaccelerated by 100× in energy, can contribute synchrotron luminosity to cluster radio halos, relics, and phoenices. Reacceleration may be aided by CRe scattering upon the ion-Larmor-scale waves that spawn when ICM is compressed, dilated, or sheared. We study CRe scattering and energy gain due to ion cyclotron (IC) waves generated by continuously driven compression in 1D fully kinetic particle-in-cell simulations. We find that pitch-angle scattering of CRe by IC waves induces energy gain via magnetic pumping. In an optimal range of IC-resonant momenta, CRe may gain up to ∼10%–30% of their initial energy in one compression/dilation cycle with magnetic field amplification ∼3–6×, assuming adiabatic decompression without further scattering and averaging over initial pitch angle. 

We study the linear stability of a planar interface separating two fluids in relative motion, focusing on the symmetric configuration where the two fluids have the same properties (density, temperature, magnetic field strength, and direction). We consider the most general case with arbitrary sound speed cs, Alfvén speed vA, and magnetic field orientation. For the instability associated with the fast mode, we find that the lower bound of unstable shear velocities is set by the requirement that the projection of the velocity on to the fluid-frame wavevector is larger than the projection of the Alfvén speed on to the same direction, i.e. shear should overcome the effect of magnetic tension. In the frame where the two fluids move in opposite directions with equal speed v, the upper bound of unstable velocities corresponds to an effective relativistic Mach number $M_{\rm re}\equiv v/v_{\rm {f}\perp }\sqrt{(1-v_{\rm {f}\perp }^2)/(1-v^2)} \cos \theta =\sqrt{2}$, where $v_{\rm {f}\perp }=[v_{\rm {A}}^2+c_{\rm s}^2(1-v_{\rm {A}}^2)]^{1/2}$ is the fast speed assuming a magnetic field perpendicular to the wavevector (here, all velocities are in units of the speed of light), and θ is the laboratory-frame angle between the flow velocity and the wavevector projection on to the shear interface. Our results have implications for shear flows in the magnetospheres of neutron stars and black holes – both for single objects and for merging binaries – where the Alfvén speed may approach the speed of light.

 

Turbulence driven by active galactic nuclei activity, cluster mergers, and galaxy motion constitutes an attractive energy source for heating the intracluster medium (ICM). How this energy dissipates into the ICM plasma remains unclear, given its low collisionality and high magnetization (precluding viscous heating by Coulomb processes). Kunz et al. proposed a viable heating mechanism based on the anisotropy of the plasma pressure under ICM conditions. The present paper builds upon that work and shows that particles can be heated by large-scale turbulent fluctuations via magnetic pumping. We study how the anisotropy evolves under a range of forcing frequencies, what waves and instabilities are generated, and demonstrate that the particle distribution function acquires a high-energy tail. For this, we perform particle-in-cell simulations where we periodically vary the mean magnetic fieldB(t). WhenB(t) grows (dwindles), a pressure anisotropyP_⊥>P_∥(P_⊥<P_∥) builds up (P_⊥andP_∥are, respectively, the pressures perpendicular and parallel toB(t)). These pressure anisotropies excite mirror (P_⊥>P_∥) and oblique firehose (P_∥>P_⊥) instabilities, which trap and scatter the particles, limiting the anisotropy, and providing a channel to heat the plasma. The efficiency of this mechanism depends on the frequency of the large-scale turbulent fluctuations and the efficiency of the scattering the instabilities provide in their nonlinear stage. We provide a simplified analytical heating model that captures the phenomenology involved. Our results show that this process can be relevant in dissipating and distributing turbulent energy at kinetic scales in the ICM.